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Abstract. We have investigated the dynamics of bound particles in multilevel current-carrying 
quantum dots. We look specifically in the regime of resonant tunnelling transport, where 
several channels are available for transport. Through a non-Markovian formalism under the 
Born approximation, we investigate the real-time evolution of the confined particles including 
transport-induced decoherence and relaxation. In the case of a coherent superposition between 
states with different particle number, we find that a Fock-space coherence may be preserved 
even in the presence of tunneling into and out of the dot. Real-time results are presented for 
various asymmetries of tunneling rates into different orbitals. 



1. Introduction 

Some of the most peculiar features of quantum theory, such as the existence of quantum 
superpositions and of entangled states, are typically destroyed by uncontrolled and ultimately 
inevitable interactions with a surrounding, often incoherent environment [I]. 

Most of the recent theoretical analysis done in the area of low-dimensional dynamical quantum 
systems has been either for open Markovian systems, where the past memory of the system is 
neglected, or for closed unitary systems, where the dynamics are reversible. Furthermore, much 
of the research done until recently has focused on steady state phenomena where a Markovian 
approach can be expected to provided reasonable results. 

However, the transient behavior of the system carries a tremendous amount of information 
in the form of coherence and relaxation dynamics. Ultrafast laser pulse excitations, for example, 
are providing insight into the heterostructure dynamics on a femtosecond timescale [2j ■ A theory 
which accounts for quantum dynamical behavior on the same time scale is therefore of great 
benefit not only to experiment, but also to basic understanding. 

The formalism of non-equilibrium Green's functions (NEGF) and formalisms based on 
NEGF, have been extraordinarily successful in the analysis of much phenomena in mesoscopic 
systems O HI [5] . Notably, work in both elastic and inelastic transport in quantum dots has been 
successfully done by means of NEGF or derivatives of this formalism [6j [7]. However, NEGF 
is inherently a closed-system formalism where the system has Hamiltonian dynamics [8J, and 
thus does not account for irreversibility arising from interactions with an unseen, unknown, or 
otherwise intractable environment. Promising attempts have been made to extend NEGF to 
open quantum system, for example, by treating the environment as a correction to the system's 
self-energy [9], or by calculating two-time correlation functions [10], effectively separating the 
time scales into transient and steady-state regimes. These and other approaches constitute 



significant progress toward a full non-Markovian open system analysis within the NEGF 
formalism. 

Even so, if the evolution of the (possibly correlated) many-body system states themselves are 
sought, one must calculate the Green's functions and then extract the density matrix, which in 
the non-equilibrium case may not be unique |10| . An alternative formulation, one we adopt in 
the sequel, is to work with the density matrix directly. This arguably provides a more intuitive 
description of the dynamics of non-Markovian open quantum systems, and it readily yields the 
actual states of the system as well as the quantum coherence between them. 

The focus of this paper is to thus present a generic formalism to investigate the dynamics of 
a few-electron quantum dot in a fully non-Markovian fashion in the density matrix formalism. 
Section [2] presents a brief outline of the derivation of the evolution equations for the density 
matrix of the system in both the Markovian and non-Markovian regimes. Section [3] presents 
the model Hamiltonian for our system as well as the analytical expression for the transition 
tensor for this model. Section [4] illustrates the non-Markovian approach by considering the 
case of a multilevel dot in a regime where two channels are available for transport. Results are 
presented for the transient population probabilities in regimes where the tunnelling strength 
between system and leads varies from orbital to orbital. In Sec. [5] we give a brief discussion on 
Fock-space coherences and their implications for Quantum Information. Finally, in Sec. [6j we 
present a summary and conclusion. 

2. Evolution of the density matrix for an open quantum system 

In isolated quantum systems, the time evolution can be described by a single state vector 
\ip(t)). Once the system is coupled to a larger environment whose degrees of freedom are either 
unknown, intractable, or uninteresting, the environmental degrees of freedom may be traced out 
and the resulting state of the system alone is described not by a state vector, but by a density 
operator [TT] p{t) = ^2 n Wni^lipn^)) (ip n (t)\, where W n (t) is the probability that state \tfj n (t)) is 
occupied at time t. The time-dependence of the density operator can be obtained by an iterative 
integration of the Liouville-von Neumann equation, which gives 

p(t) = -C(t)p{0) + f dt'£(t)C{t')p{t') + fdt' f dt"C(t)C(t')C(t")p(t") + ..., (1) 
Jo Jo Jo 

where C(t) is a Liouville superoperator such that, C(t)F = i?i _1 [if(i), F] for any operator F 
and Hamiltonian H(t). 

Equation ([TJ represents the evolution of the complete system plus environment. Under the 
Born approximation, the evolution equations are truncated at second order in the interaction 
Hamiltonian. After some manipulation the resulting equations for the matrix elements of 
the density operator read, 

Pab(t) = V f dt' Pcd {t')R abcd {t - 0e iW ', (2) 

c,d J ° 

where, R a bcd(t — f) is a Redfield-type transition tensor containing all the information 
characterizing the system, the environment, and the coupling between them, and where 
labcd = w„(, — uj c d = (E a — Ef, + E^ — E c )/h denotes energy differences between system states. 

For an environment with a much greater number of degrees of freedom relative to those of 
the system, and for a weak enough system-environment coupling, we make the assumption that 
the system has no effect (no back-action) on the environment, that the system is not correlated 
with the environment, and that the environment is and remains in equilibrium [11 . In this case, 
the total density matrix for the system plus environment can be written as, 



p{t) Kt p S ystem(t)0R(O). 



(3) 



where p S ystem(*) is the density matrix for the open quantum system, and £?r is the density matrix 
for the environment. 

2.1. Markovian Approach 

The Markov approximation assumes an environment which relaxes on a time scale time much 
faster than the characteristic time of the system. The environment correlation functions vanish at 
such a fast rate, that the reversibility of the system is essentially destroyed, thus also destroying 
the memory of the system [I] . Thus, p S ystem(t) depends only on the present value of p{t) . That is, 
it becomes local in time. For time intervals t — t' much greater than the environment's correlation 
time r, the correlation functions for environment operators rapidly become uncorrelated and 
decay to zero, {F ] k {t - t')F k ,) « {F ] k {t - t'))(F k/ ) « 0. In the limit t - t' » t the upper 
integration limit in ^ can be extended to infinity with negligible error in the calculations. 
Since the time when the coupling is turned on is arbitrary, the lower integration limit can be 
taken to negative infinity, and the coupled set of integro-differential evolution equations become 
a coupled set of first order ordinary differential equations: 

/oo 
dt'R abcd (t - tV W ' = Yl Pcd(t)W abcd . (4) 

ca -°° cd 

This in turn leads to a Fermi's Golden rule for transitions with strict energy conservation for all 
transitions. 

Although much of the recent theory in the area of non-equilibrium quantum dot systems 
has made use of Markovian-type approximations p2J [131 El El HHJ E], the approach has 
fairly severe consequences on the transient behavior of the system, especially when considering 
coherence properties of the system at short time scales. Nonetheless, the Markov approximation 
may be expected to yield reasonable results at sufficiently long times [TT] , when the system has 
reached a steady state, or when considering intermediate time scales of averaged behavior, such 
as the average total current through a system. 

2.2. Non-Markovian Approach 

In order to account for the transient dynamics, especially for the characteristics of the coherence 
between system states, memory effects must be preserved. In our approach we make use of the 
convolution theorem for Laplace transforms, to represent Eq. ^ as a set of algebraic equations 
in Laplace space, 

SPab{s + i^ab) ~ p(0) = ^ Pcd ( S + iuJ cd)Rabcd{s + lU) ab ). (5) 

cd 

This set of coupled algebraic equations can be solved analytically for p a b(s) for only a few 
system states (~ 2-4), since the computational effort rapidly increases with the number of 
available transport channels. In the regime of sequential resonant tunnelling transport through 
a quantum dot containing non-interacting electrons, the transport channels are those single- 
particle system states \a), whose energy lies within the bias window. These single-particle 
states define the possible dynamical many-body states of the system as those involved when 
the single particle states are empty or occupied. Thus, for k channels, a minimum of 2 k many- 
body system states are required (for empty or occupied). Furthermore, (2 fc ) 2 density-matrix 
elements are required do describe the population probabilities as well as the coherence between 
the states. The temporal evolution of the density-matrix elements themselves are governed by 
(2 fc ) 4 transition tensor elements (Redfield-type in the Born approximation). For larger sizes, 
a numerical approach may be utilised to invert the linear system in Eq. ^ at each (complex- 
valued) point s, thus obtaining p ab (s). The solutions p a b(s) in Laplace space can then be 



brought back to time space by applying an inverse Laplace transform. This is done numerically 
by performing a Bromwich integral [TH] over a suitably chosen contour. For every point s along 
the Bromwich contour, a matrix inversion is performed to evaluate all components of p a b{ s )- 
Such a Bromwich integration is (numerically) performed [19] at each time step, thus obtaining 
Pdb(t). 

2.3. Violation of Positivity 

It has been established by Kraus [20], Lindblad [21] . and Gorini, Kossakowki, and Sudarshan [22], 
that the evolution of an open quantum system must remain completely positive to ensure the 
physical validity of the states of the system at all times. However, this positivity may be broken 
by introduction of ad hoc relaxation rates, by the choice of basis states, or by the truncation of 
higher-order interaction terms in the evolution equations (such as in the Born approximation), 
which can yield serious inconsistencies such as negative probabilities. Despite this, much of the 
investigations of non-Markovian dynamics acknowledge that for a set of initial conditions, the 
theory will reproduce physically consistent results. Our findings agree with this premise [23 . 

3. Model 

We model our system as a single generic quantum dot coupled to a source and a drain reservoir 
by a tunneling Hamiltonian [24 . The total Hamiltonian is given by, 

H = H S + H QD + H D + H T , (6a) 

where H$ and Hp are the source and drain Hamiltonians, respectively. These are taken to be 
non-interacting Fermion systems shifted by the bias: 

H S(D) = J2^ s ( d ) ± \ eVB } d \{d) d s{d)i ( 6b ) 

s{d) 

with dK* a creation operator for the source (drain), and d s ^ an annihilation operator. 
The Hamiltonian for the quantum dot in Eq.(|6a[) is given by, 

Hqd = y~](fc^ + eV g )c\ci + Vint, (6c) 

where the single-particle energies hcoi are all shifted by the applied gate voltage, and Vi n t is the 
interaction among the confined particles. 

The tunneling Hamiltonian describing the reservoir-dot coupling is given by, 

H T = Yl { T k d l c k + h.c^ , 

k,r=(s,d) 

where T?^ is an energy-independent tunneling coefficient for a particle tunneling from the 
single-particle state \k) in the dot to the source (drain) reservoir. Finally, we assume zero 
temperature such that the reservoirs are in their respective non-degenerate ground states. 

As mentioned above, all system and environment information, including the characteristics of 
the interaction between them are contained in the memory kernel, R a bcd{T~) appearing in Eq.Q. 
For the model in Eq. (|6]), we obtain 

+ Kfar) {{% Bi » R TKL - i^TKL} , (7) 



(6d) 



where the indices a, b, c, and d denote many-body states, a, and f3 are denote single-particle 
states, and R = S, D runs over the source and drain leads. The energy 4>t(b) denotes the top 
(bottom) of the band, and fj, R is the chemical potential of reservoir R. Furthermore, we have 
defined 

= e l ^ aT - e iuJ ^ T , (8a) 
Kfcd = (a\ct\c)(d\c p \b) - (a\c)(b\cj p \d), (8b) 

and 

<^) = ^(T*Te) R , (8c) 

where u x>a = (E x — E a )/Ti is a frequency denoting the energy difference between states \x) and 
| a), T a is the tunneling coefficient to a single particle state \ot) in the dot, Nr is the density of 
states for the 2DEG environment, and r is a time parameter. 

With this general form for the transition tensor, given a set of basis states, the individual 
elements of the tensor are calculated. The form of the Redfield transition tensor allows for 
independent tuning and analysis of barrier width symmetries, and of asymmetries between core- 
states and edge-states couplings to the reservoirs, by means of the independent tunneling rates. 

4. Two-channel (four-state) system 

As an illustrative example, we present results for the diagonal elements of p{t) in the above 
model. For definiteness, we consider a quantum dot with N confined particles, and with two 
resonant tunnelling transport channels available within the bias window. Each channel involves 
a particle- number fluctuation between N and N ± 1, and can involve either the ground or 
first excited state of the iV-particle system. In general, the availability of k transport channels 
involves a minimum of 2 k states. In the present case, the four relevant states are 

\0) = \(N-l) g . s .), |l) = |iV 9 . s .), \2) = \{N+l) g . s ), |3) = |(iV)„). (9) 

Here, the state |0) denotes the (N — l)-particle ground state of the system, |1) and |2) denote 
ground-states of the TV and (N + 1) -particle system respectively, and |3) denotes the first excited 
state of the iV-particle system. A schematic representation of the system is shown in Fig. [TJ 

Each orbital will in general be coupled differently to the leads, owing to the detailed shape 
of the wave function [25]. We identify in Fig. [T] edge and core orbitals. Core orbitals are weakly 
coupled to the leads owing to a poor (s-type) overlap with the lead states whereas edge orbitals 
(p-type) are more strongly coupled to the leads. 

Figure [2] shows the real-time evolution of the population probabilities of the system p n n(t) 
in the large bandwidth limit. Results are shown for a bias of Vbias = 6 meV, symmetric about 
the Fermi energy eFermi- We assume two channels within the bias window, E a = eFermi + 1 meV 
and Ep = eFermi — 1 meV. Finally, we keep the barriers fully symmetric (i.e. T' source = 2 ndram ) j 
while the ratio between tunneling to core states and tunneling to edge states, is varied (i.e. 
T a /T@ > 1). Several features are evident in Fig. [2] and we discuss three in particular. 

First, we see that for the smallest couplings to the core orbitals (the two plots on the right 
of Fig. [2]), the four levels couple into two distinct pairs. States |1) and |2) — the N and (N + 1)- 
particle ground states — are strongly coupled, as are states |0) and |3) — the (N — l)-particle 
ground state and the iV-particle first-excited state. We can understand why this occurs with 
recourse to Fig. [T] The transition between states |1) and |2) is through an edge orbital and this 
coupling is stronger than the transition between |1) and |0) which involves a core orbital. Similar 



arguments apply for the coupling between states |3) and |0) (strong coupling) and between |3) 
and 1 2) (weak coupling). 

Second, in the steady state {t — > oo), all occupation probabilities tend to the same value 
of 1/4. This is seen regardless of the tunnelling strengths of core and edge states. The equal 
probability of 1/4 for each level can be understood as being due to the symmetric barriers 
between the dot and the source on the one hand, and between the dot and the drain on the 
other. Since these barriers are symmetric, any level will have an equal probability of being either 
occupied or unoccupied. 

Third, although all four states are equally occupied in the steady state, the time to actually 
reach the steady states does depend on the relative couplings to the core and edge orbitals. The 
fact that the time taken to reach the steady state increases as the tunnelling to the core state 
decreases can be understood as the effect of decreasing the available tunnelling pathways. Fewer 
pathways available means that the system takes longer to reach the steady state. In the limit of 
zero tunnelling to the core state, for example, states |3) and |0) will never become occupied. In 
this case, the remaining two levels will each reach an occupation probability of 1/2 rather than 
1/4. 

5. Fock- Space Coherence 

Finally, we comment on the various forms of coherence available to systems like the present 
one in which the particle number fluctuates. By Hilbert-space coherence, we denote coherence 
between states with identical particle numbers — states |1) and |3) in Eq. for example. In 
contrast, a Fock-space coherence can exist between states of differing particle number. In the 
system considered above, we can have a Fock-space coherence between states differing by one 
particle — four possibilities for the states in Eq. Q — or by two particles — states |0) and |2). In 
general, the interaction and interplay between these distinct classes of coherence is governed by 
the Redfield tensor R a bcdif)- 

For the k = 2 channel system, k 2 = 4 four states are required, as outlined above. This 
yields (k 2 ) 2 = 16 density matrix elements, and thus, (fc 2 ) 4 = 256 Redfield tensor components, 
many of which are not independent and many of which vanish. In Fig. [3] we show a graphical 
representation of all non-zero components of the Redfield tensor, for the same 2-channel case 
described above for the model Eq. ([6]). In this matrix representation, we see three distinct and 
decoupled submatrices. The top submatrix represents the transitions between the population 
probabilities and all Hilbert-space coherence terms. The remaining two block describe the 
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Figure 1. Schematic rep- 
resentations of system un- 
der consideration. Four sys- 
tem states are considered 
depending on the relative 
occupation of the single- 
particle orbitals \a) and |/3), 
each of which can in prin- 
ciple have different tunnel 
couplings to the leads. 
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Figure 2. Time evolution of population probabilities in a quantum dot with 2 transport channels 
and four states. (See Eq. ([9]).) The plots are for symmetric source and drain tunnel barriers, 
and varying core edge asymmetry. We assume a 6 meV bias symmetric about the Fermi energy, 
and two transport channels at energies epsiloriF eTm i ± 1 meV. 



Fock-space coherence between states differing by either one (middle submatrix) or two (bottom 
submatrix) particles. 

The remarkable structure evident in Fig. [3] implies that each block undergoes its own 
independent time evolution. In turn, this suggests that, should a Fock-space coherence be 
established in the system at any time, this coherence can be robust and long-lived even in the 
rather severe perturbation of transport carriers being continually injected and removed from the 
system. This may have very promising consequences for quantum information, where strong, 
accessible, and long-lived coherence is an essential prerequisite for information processing beyond 
the classical limit. 

6. Conclusion 

We have developed a generic non-Markovian formalism for the real-time evolution of the density 
matrix of a quantum dot weakly coupled to source and drain reservoirs. Within the tunneling 
Hamiltonian approach, we have analytically derived an explicit transition (Redfield) tensor 
describing the full dynamics of the system. In the case of two transport channels (four states) 
we observed the real-time evolution of the population probabilities, taking account differing 
tunnelling rates to different orbitals. Two distinct types of coherence became evident from 
the results, a Hilbert space coherence between states with same particle number, and a Fock- 
space coherence between states with different particle numbers. Under appropriate conditions, 
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Figure 3. Representation of the non-zero Redfield tensor elements for the non-interacting, non- 
Markovian theory for the 2 transport-channel case with the four-states: |0) = | (-ZV — l) g . s .), |1) = 
\Ng. s .}, |2) = |(iV+ and |3) = |(A r ) e . s .). Three distinct decoupled submatrices become 

apparent: The top submatrix represents the transitions between the population probabilities and 
the coupled coherence between them (Hilbert-space coherence) , the middle submatrix represents 
transitions relating states with particle numbers differing by 1, and the bottom submatrix (a 1 
by 1 matrix) represents a transition relating states with particle numbers differing by 2. 



the Fock-space coherence is decoupled from the evolution of the rest of the system, suggesting 
promising consequences for quantum information processing. 
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